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, The top quark pair production and decay are considered in the framework of the smeared- 

; 1 mass unstable particles model. The results for total and differential cross sections in 

' vicinity of tt threshold are in good agreement with the previous ones in the literature. 

The strategy of calculations of the higher order corrections in the framework of the 
model is discussed. Suggested approach significantly simplifies calculations compared to 
the standard perturbative one and can serve as a convenient tool for fast and precise 
preliminary analysis of processes involving intermediate time-like top quark exchanges 
in the near-threshold region. 
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1. Introduction 

The top pair production and decay are the key processes for precision tests of the 
Standard Model (SM) (see e.g. Ref. [1] and references therein). They were intensively 
studied in the framework of the Quantum Chromodynamics (QCD) and Electro- 
Weak (EW) perturbation theory during last two decades, and various methods and 
schemes were proposed. The major goal of these investigations is to define the basic 
physical parameters of the top quark, such as its mass, width and couplings with 
other SM particles. In the past, the top quark physics was one of the primary 
research objectives at Tevatron. Nowadays, the biggest attention is paid to the 
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process of the top quark production at the LHC (see e.g. Refs. However, the 

highest precision measurements of the top quark properties can best be reached at 
the future Linear Colhder (LC) which supposedly operates in a clean experimental 
environment. The top quark physics is one of the most interesting and chahenging 
targets for future e+e~ or /x^/x" LC experiments [4j[5]. 

The top pair production is followed by a decay chain with intermediate gauge bo- 
son states, i.e. the full process under consideration is e+e~ — > t*t* hbW^W~ — >■ 
hhAf . The widths of both the top quark and the M^-boson are large, and one nec- 
essarily needs to take into account corresponding Finite- Width Effects (FWE). In 
the framework of the standard perturbative approach, these effects are typically 
described by means of dressed propagators which are regularized by the total decay 
width. In order to analyze the full process of the top pair production relevant for 
phenomenological studies, we also have to take into account the background con- 
tribution coming from many other topologically different diagrams leading to the 
same six-fermion final states, which is a rather non-trivial task. 

The Born-level cross-sections of the processes e^e" — > bbudii^Vf^ and e^e~ 
bMq were calculated in Refs. Oil] and [8], respectively. Other exclusive reactions 
with bbdu^^Vf^i, bbcsdu and bb^'^v^T~Vr final states were considered in Ref. [5]. In 
particular, it was shown that the contribution of the top-pair signal e+e" — >■ t*t* — 
664/ is dominant, but the background (caused by one- resonant or non-resonant 
diagrams) can be quite significant too. However, it can be drastically decreased by 
applying certain kinematical cuts on the appropriate invariant masses. 

The QCD corrections for the reaction e^e~ — ^ tt in the continuum above 
the threshold were previously obtained in Refs. [HI HD]- As well as the one- loop 
EW corrections were calculated in many papers (for corresponding references, see 
e.g. Introduction in Ref. [H]). Concerning radiative corrections (RC) to reaction 
e~''e~ — >■ 664/ with six-fermion final states, the situation is more complicated and 
less clear ill . At the tree level, any of the reactions receives contributions from sev- 
eral hundreds of diagrams. The calculations of the full 0{a) radiative corrections 
are very complicated, and different approximation schemes are typically applied. 
The most detailed analysis of the exclusive reactions e^e~ — ?> bb^^v^^^^D^ and 
e^e~ — ^ bbdu^~D^ was performed in Ref. [llj . In this paper, the cross-sections were 
calculated taking into account the leading radiative corrections, such as the initial 
state radiation (ISR) and factorizable EW corrections to the on-shell top-pair pro- 
duction, to the decay of the top quark into bW and to the subsequent decays of the 
M^-bosons. Usually, such calculations are carried out automatically by Monte Carlo 
techniques (see Ref. [11] and references therein). 

In this work, we consider reactions like e+e^ -> t*t* — ?> 664/ with any four- 
fcrmion final states 4/. The analysis is performed in the framework of the smeared 
mass unstable particles model (below, SMUP model) [121 IE] ■ Due to exact factor- 
ization at intermediate i, i and W'^ , W~ states, the cross-section can be represented 
in a simple analytical form which is convenient for analytical and numerical anal- 
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Fig. 1. Feynman diagram of the top quark signal process e+e — > ft* — > bbW~^W —> bbAf. 

ysis. So far, we have appHed the SMUP approach only for unstable gauge boson 
productfon and decay (see e.g. Refs. [TH [TSl [H] ) . As a continuation of our earlier 
studies, in this work we test the SMUP approach for the case of unstable fermions, 
i.e., specificaUy, top quarks. In our calculations, we take into account NLO radiative 
EW and QCD factorizable corrections which dominate close to tt threshold. Also, 
we illustrate the influence of the mass smearing effects and various radiative correc- 
tions (RC's) on the differential cross-sections. The results are compared with ones 
calculated by using the standard perturbative methods pj], where cross-sections 
were represented for case of full 2 — )■ 6 process and, separately, for the top signal 
contribution alone. It was shown that in the Born approximation the results coin- 
cide with a rather high precision, and deviations of the higher-order corrected results 
from the standard ones are at the percentage level. So, the suggested approach can 
be applied in a fast preliminary analysis of various complicated processes involving 
intermediate top quark exchanges in the Standard Model and beyond. 

Note, here we do not consider the near-threshold effects caused by the generation 
of the coupled ti state, which were considered in detail in many previous studies 
(see, for instance, Ref. [17] and references therein). We postpone this issue for a 
forthcoming study. 



2. The model cross-section of the top-pair production and decay 
at the tree level 

The process of top-pair production with subsequent decay e+e^ — >■ t*t* — > 
bbW^W^ —T' bMf is schematically represented in Fig. [TJ The full process contains 
two steps with unstable intermediate time- like states, namely, t^t and W^,W~ 
states. In this case, as was shown in Ref. [13], the double factorization takes 
place and can be described in the framework of the SMUP model [12]. Due to 
this factorisation, the full process can be divided into three stages: e+e~ — >■ t*t*, 
t*t* bbW'^W~ and W+W~ 4/. Here, the top-quarks and W-bosons are 
treated as unstable particles, and finite- width effects should be taken into account. 
The SMUP model cross-section of the first reaction e+e^ t*t* can be written 
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as 



(T(e^e^ — > = / / a{e^ ^ t{mi)t{m2))pt{i'ni)pt{rn2)dm\dm\, 

(1) 

where mo ~ 2Mh [Mb is the bottom quark mass) is the threshold value of the top 
mass variable, tT(e+e^ t[mi)t{m2)) is the cross-section of top pair production 
with random masses mi and m2 and pt{m) is the probability density which describes 
the mass smearing of top quarks. In our calculations we take it in the Lorentzian 
form as [T2] 

( X 1 mrt(m) 

m = -— ^ (2) 

TT (m^ — Mj )^ + m^I f (m) 

where Tt{m) is the total decay width of the top quark with mass m. The decay mode 
t — >■ bW has a very large branching ratio Br(i — > 6W^) « 0.999, so formula ([T]) almost 
exactly describes the cascade process e+e^ — > t*t* — ?> bbW~^W^ in the stable VF- 
boson approximation. In order to take into account the instability of VF-bosons we 
have to express the top quark width r((m) « r(t — )■ bW) in Eq. as a function of 
smeared VK-boson mass r{t — ^ bW{M\Y)) with averaging over Mw- Thus, the model 
cross-section of the full inclusive process e+e^ — >■ t*t* — )■ bbW^W" bb'Y^^Af 
depicted in Fig. [T] has the following convolution form: 

/•s <.(v^-mi)^ 

aie-^e^ ^ bb^if) = / a{e+ ^ t{mi)t{m2))x 

pt{mi,m-w+)pw{mw+)dm^+x (3) 

(mo-Afb)2 

Pt {m2 ,mw-)pw {mw- )dm^w- dmldml , 

where pw(jm) is defined by Eq. ([2]). In order to describe an exclusive reaction 
e+e" — >■ t*i* — > 66/1/2/3/4 we have to replace the total decay widths of M^-bosons, 
which enter the numerator in Eq. ([5]), by corresponding exclusive ones (see Section 
4). The same result can be obtained exactly if one calculates the cross-section of this 
process explicitly in the framework of the SMUP model by using dressed propaga- 
tors of unstable particles (UP's). In Ref. [T^ it was shown that exact factorization 
of a decay chain process with UP's in an intermediate state takes place when we 
exploit the model effective propagators for fermion and vector UP's in the following 
form 

Dil) = ^J^y D,Aq)--^ p^^^^ , (4) 

where Ppiq) and Pv{q) are the denominators of the fermion and vector boson 
dressed propagators, which contain corresponding total decay widths. The structure 
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Fig. 2. The cross-sections of the processes e+e — > tt, e+e — > bbW~^W and e+e — > bb'^j, if. 

of numerators in Eq. (j4]) provides exact factorization and leads to a convolution- 
like expression ([3|) for the cross-section. So, there is a self-consistency between the 
model and UP effective theory description of the processes with UP in intermediate 
states. Thus, the process with a six-particles final state shown in Fig. [T] is described 
by a simple analytical expression (jS]) with four integrations over smeared unstable 
top and VF-boson masses. Note, that the standard perturbative treatment of the 
six-particle final states in general case leads to7V = 3-6— 4=14 independent param- 
eters, from which 13 parameters have to be integrated over [15]. Such a complicated 
problem can be solved by using involved Monte Carlo numerical simulations only. 

The results of the SMUP model calculations are presented in Fig. [21 Here, the 
dotted line represents the cross-section of the top-pair production in the stable 
particle approximation (SPA), i.e. without smearing of the top mass. The dashed line 
is the cross-section incorporating the top mass smearing or top finite-width effects 
(FWE) only, and the solid line gives the full mass smearing result, both top-quarks 
and VF-bosons. Note, that the second case corresponds to the standard treatment 
of the process e^e~ — )■ t*i* — bbW^W~ in the stable VF-boson approximation, and 
the third case - to the full process shown in Fig. [1] 

From Fig. [21 one can see that the contribution of the top quarks' FWE's is signif- 
icant (up to a few percents in the near-threshold region), while the contribution of 
VF-bosons' FWEs is small. The comparison of our results with ones in the standard 
perturbative treatment shows that deviations are typically very small. For instance, 
it was obtained in Ref. [19], that a{e+e- t*t* bbW+W^) for = 500 GeV 
is equal to 629 fb for Mt = 150 GeV and 553 fb for Mt = 180 GeV. For the same 
input data, we have obtained 630 fb and 554 fb, respectively, which are in a good 
agreement with the result mentioned above. This comparison proves the applicabil- 
ity of the SMUP model fermion propagator given by the first expression in Eq. Q . 
In Section 4, we make such a comparison for exclusive processes as well where both 
SMUP model fermion and boson propagators Eq. (|4|) are used. 

It should be noticed also that we consider the FWE's, which are significant in 
the near-threshold region, but we do not include near-threshold effects caused by 
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possible intermediate ti bound states. Since the top mean lifetime is considerably 
shorter than the hadronisation time, the bound state effect has no sharp resonant 
nature. However, it can be comparable with FWE's or mass-smearing effects under 
consideration, and this problem will be considered in more detail elsewhere. 

3. Factorizable corrections to the cross-section 

As it was shown in previous papers [11], the EW and QCD corrections give large 
contributions to the cross-section of the top-pair production at energy scales close 
to its threshold. In this Section, we describe the strategy of our model calcula- 
tions and give the total cross-section including the principal part of NLO EW and 
QCD corrections. Note, that the strategy of calculations and the choice of input 
parameters are mainly caused and defined by the effective character of the model 
treatment. In the framework of the SMUP model, the instability (or finite width) of 
unstable particles is accounted for by the smearing of their masses, i.e. by the prob- 
ability density function p{m). In turn, this function contains momentum dependent 
parameters M(q) and T{q) in analogy with the standard perturbative treatment 
which uses dressed propagators. So, in that sense the corrections of self-energy type 
are already included at the "effective" tree level, and it is reasonable to use an ef- 
fective couplings, such as running coupling, absorbing the major part of vertex-type 
corrections. 

In our calculations we have used the following input data I20j : 

a{Mz) = 0.00781763, as{Mz) = 0.118, sin^ 9w{Mz) = s| = 0.2313, 

Mz = 91.1876 GeV, Mw = 80.399 GeV, Mt = 172.9 GeV. (5) 

The running coupling constants afc(Q^), fc = 1,2,3 were used in the one-loop ap- 
proximation: 

- l-(/3.7t)Tn(QVM|) ' - 
The cross-sections are calculated including the following corrections: 

• Vertex and self-energy type corrections for stable particles are mainly in- 
cluded into running couplings 

• Self-energy corrections for unstable particles are included into the proba- 
bility density function p{m), which describes the smearing of UP's masses. 

• Initial state radiation (ISR) is described by the photon radiation spectrum 
pT| [22] , and the bremsstrahlung from the final t-quark states - by vertex 
Q-dependent factor [23]. 

• QCD corrections to the top production and decay are described by the ver- 
tex multiplicative factor [23] . 
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Fig. 3. The higher order corrected cross-sections of the process e+e — >■ bb'^j: if. 

• Contribution of the box diagrams to the total cross section was evaluated 
at energy scales close to the threshold by using numerical FormCalc v7.3 
[24] routines. 

The higher order corrected cross-sections of the inclusive process e+e" — )■ t*t* — 
hh'Y^jAf are shown in Fig. |3l There, the dotted line represents the Born model 
cross-section, the dashed line - the cross-section with ISR and the solid line - the 
cross-section with total factorisable corrections (without box diagrams contribu- 
tion). From the figure, one can see that the main contribution is given by ISR 
correction, which significantly reduces the cross-section in the near-threshold en- 
ergy range and increases it at energy scales above ~0.6 TeV. At large energies 
(y^ > 0.5 TeV) the contribution of EW and QCD corrections becomes significant 
and has to be properly taken into account. 

In Fig. |4] we present the invariant mass distribution and illustrate the influence 
of various corrections on it. One can see that the corrections which we have taken 
into account according to the procedure above give noticeable contribution into 
this distribution in the peak area. We, also, illustrate such influence on the angular 
differential cross-section presented in Fig. [S] Again, we notice that this influence is 
quite significant and should be taken into account. 

4. The cross-sections of exclusive processes 

So far, we have considered the cross-section of inclusive process e"'"e~ — t*i* — >■ 
bbJ2f^f where the final state is summed up over all possible fermion fiavors. As 
was noticed in the second Section, in order to get the cross-section of exclusive 
process e+e^ t*i* — > fo6/i/2/3/4 we can include the corresponding branching 
ratios Br(Vl^ -> /1/2) and Br(Vl^ -> fzli)- Acting this way we obtain 

a(e+e- ^ bbhfihh) = o^e^e- ^ 66^4/)Br(I^ ^ /i/2)Br(I^ ^ /aA) . (7) 

/ 
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Fig. 4. Invariant top mass distribution. 



where we omit intermediate virtual t*t* state for simplicity. This relation directly 
follows from the Eq. ([3|) when one substitutes a partial decay width of the W- 
boson into numerator of the probability distribution function pw {m) instead of the 
total width. It can be also derived by straightforward calculation of the a{e^e'^ — 
bhfif^fzfi) in the framework of the effective theory (see Ref. [13^). 




-1.0 -0.5 0.0 0.5 1.0 



Fig. 5. The angular differential cross-sections for the process e^e — > bb J^j- 4/. 



The expressions for the branchings ratios Br(T/F — !■ /1/2) were considered in 
detail in Ref. i23:. Here, we use very simple but sufficiently precise formulae which 
incorporate QCD corrections: 

^^(^ '^^^ - 9(l + 2a.(M.)/3.) ' ^'^(^ ^ = 3(l + 2a.(M,)/3.) ' 

(8) 
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where Vik are elements of the Cabibbo-Kobayashi-Maskawa mixing matrix. We, 
also, employ the QCD corrected expression for the top quark width [23l [25]: 

Tit ^ bW) = l^a^iMtWtbl^ VQCD Mt f{Mu Mw,Mb), (9) 



where 



2A2\ 1/2 



2as{Mt) /27r2 5 



??QCD = 1 - 



Stt V 3 2 



Using Eqs. ([7|)-(l9l) we can calculate the exclusive cross-section for an arbitrary 
six-fermion final state (6^/1/2 /3/4)- Such calculations taking into account the factor- 
izable EW corrections were performed within the standard perturbative approach 
for the case of {bbfi^Vf_,fi~D^) and final states in Ref. [H]. In this work, 

the full set of topologically different Born diagrams leading to the same six-fermion 
final state was considered. It was shown, that certain cuts on invariant masses of 
the Wb and fifk pairs, which correspond to intermediate t,t and states 
for signal diagrams, significantly reduce the relative contribution of the background 
(see Table [J). 

In Tables [U and [2] the cross-sections are given for two distinct reactions 

(1) : e^e^ bb^'^v^ji'v^, (2) : e^e^ — > bb^'^Vi^idu. (11) 

for the energies -/s — 430, 500, 1000, GeV. In Table [1] the cross-sections are pre- 
sented in the Born approximation for total set of diagrams (c^^j.^ (total)) and for the 
signal diagrams (c^ornC**^*))' where fc = 1, 2 denotes the first and second reactions 
in Eq. (ITT|) . respectively. These values (in fb) are taken from Table 1 in Ref. [TT] 
and are calculated with the kinematical cuts 5i < 0.1, where 5i is the deviation of 
the ratio m™''" /Mi from unity and index i is related to different t, t, W~ states 
(for more details, see Ref. |11)). 



Table 1. Born-level cross-sections of the processes (1) and (2) in Eq. (Illl l. 



yi, GeV 


•^si-n (total) 


o-r' (t*t*) 


(total) 


f^R^ (t*f*) 

Born^ ' 


430 


5.9117 


5.8642 


17.727 


17.592 


500 


5.3094 


5.2849 


15.950 


15.855 


1000 


1.6387 


1.6369 


4.9134 


4.9106 



In Table [2] the results for the total cross sections (in fb) of processes (1) and 
(2) from Eq. (|lip in the Born approximation are shown in the second column. 
The cross-sections with separate ISR and factorizable EW (FEWC) corrections are 
presented in the third and forth columns, respectively, and the cross-section with 
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both the FEWC and ISR corrections included - in the fifth column. All values are 
calculated with the kinematical cuts mentioned above. 



Table 2. Comparison of the exclusive cross-sections of Ref. [11] and ones obtained in the present 
work. 



V3, GeV 


"Born 


0"Born+ISR 


""Born+PEWC 


0"Born + ISR+FEWC 


e+e — !> bUfi/i'^bfi i/fj,, Ref. 11 


430 


5.8642(45) 


5.2919(91) 


5.6884(55) 


5.0978(53) 


500 


5.2849(43) 


5.0997(51) 


4.9909(49) 


4.8085(48) 


1000 


1.6369(15) 


1.8320(18) 


1.4243(14) 


1.6110(16) 


e+e — > bu^fi^bfj, u^, this work 


430 


5.86476 


5.27613 


5.77727 


5.19941 


500 


5.27352 


5.08651 


5.18407 


5.00291 


1000 


1.63061 


1.83508 


1.58925 


1.79079 






e+e^ — fei/ji/x+ 


bdu, Ref. 11 




430 


17.592(13) 


15.857(20) 


17.052(16) 


15.283(16) 


500 


15.855(13) 


15.311(15) 


14.977(16) 


14.438(14) 


1000 


4.9106(46) 


5.4949(55) 


4.2697(40) 


4.8287(47) 


e+e — > buij,ii~^bdu, this work 


430 


17.8163 


16.0351 


17.5540 


15.8019 


500 


16.0203 


15.4517 


15.7516 


15.1979 


1000 


4.95397 


5.57465 


4.82889 


5.44011 



From Table [H it follows that the differences of the model and standard Born 
cross-sections are of an order of 0.1 percent and ISR correction increases it only 
slightly. In principle, these deviations can be further reduced. The situation becomes 
worse, when we take into account all major corrections. The deviations increase 
and become up to a few percents. This discrepancy is caused by the fact that in 
Ref. |11| an additional contribution from the non-signal (background) diagrams 
was included while we consider the signal contribution only. Moreover, we do not 
include the contribution of the box diagrams which becomes very important at 
large energies far from the threshold. According to estimations in the framework 
of the standard perturbative treatment, the box diagrams contribution is of an 
order of a few percents in the near-threshold energy range. Rough estimations in 
the framework of the SMUP model give the box contribution equal to 1.5 — 2 
percents in the energy region under consideration, and these estimations decrease 
the deviations. However, in the framework of the SMUP model, as well as in the 
effective theory of UP, the higher order corrections have an effective character, and 
a consistent formulation of the perturbative treatment with this model is required. 
This problem leads to using the model propagators inside loop diagrams, but the 
validity of such a procedure has still to be justified theoretically. In particular, one 
should first analyze the asymptotic properties of the propagators. The analysis is 
not carried out yet, but it is in progress now. Note, that a good agreement of the 
SMUP model and standard Born-level results, illustrated in Table [21 provides a 
good basis for such an analysis. 
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5. Conclusion 

The production of the ti pair and its subsequent decay into six fermion final states 
in e^e~ annihilation has been previously analyzed within the standard perturbative 
treatment in a vast literature. In this work, we performed the corresponding analysis 
in the framework of SMUP model. So far, this approach was applied mainly to the 
gauge boson production, where the structure of the model boson propagators was 
successfully tested [TH [151 [H] . In the present work, we have tested the structure 
of the model fermion propagator, and the top quark production mechanism has 
been chosen as an important example. It was shown that the results of Born-level 
calculations are in a good agreement with the standard perturbative ones, providing 
the applicability of the SMUP approach to the top-quark production and decay 
processes. 

The SMUP model provides simple analytical expressions for the total cross- 
sections of inclusive and exclusive processes with top quark pair production and its 
subsequent decay. It is a convenient and simple instrument for description of com- 
plicated multi-step processes with unstable particles participation. The precision of 
this approach at the tree level is of an order of 0.1 percent or better. The method 
gives a possibility to include, in principle, all factorizable corrections. Our approach 
can be useful in a preliminary analysis of complicated processes with intermediate 
time-like top quark exchanges within the Standard Model and beyond. 
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